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Clean predictions are presented for all the spin-averaged heavy-light meson spectroscopies. A new 
symmetry is identified wherein the energy eigenstates have a universal dependence on both the light 
and heavy quark masses. This universality is used in an efficient analysis of these mesons within 
the QCD string/flux tube picture. Unique predictions for all the D, D s , B, and B s type mesons in 
terms of just four measured quantities. 



I. INTRODUCTION 

One of the most promising and least developed areas of 
hadron spectroscopy is the excited heavy-light (HL) meson. 
Although at present only a few states of each flavor have 
been observed future discoveries at B factories, CLEOc, 
HERA, and hadron colliders will surely change this situ- 
ation. In addition to the ordinary qQ states we should 
observe hybrid and possibly multi-quark confined mesons. 
It is important therefore to reliably predict where the stan- 
dard HL mesons lie and to explore the close relationship 
, between the D,D S ,B and B s families of HL mesons. 

A striking observed fact for HL systems is that hyper- 
fine splittings are independent of light quark flavor. For 
example 

D* s -D s ~ D*-D ~ 142 MeV, (1) 
B* s - B s ~ B* - B ~ 46 MeV . (2) 

This apparent lack of light quark mass dependence in these 
differences is certainly not that expected in the popu- 
lar Breit-Fermi type semi-relativistic interaction which is 
proportional to the inverse product of the quark masses 
(1/toqto). In the following we show that this is an ex- 
ample of a larger "universal light quark mass dependence" 
(UMD) ultimately a consequence of relativistic kinematics. 
We note in passing that the ratio of the above hyperfine dif- 
ferences docs however reflect the inverse ratio of the heavy 
quark masses. 

We start by proposing and supporting the concept that 
all HL energy eigenstates have the same light quark mass 
dependence and hence all differences containing the same 
light flavor are independent of light quark mass. We take 



this as an organizing principle to analyze the various HL 
systems. In particular we find a functional relation be- 
tween strange and non-strange light quark masses, m s and 
m Ui d, and we predict the spectra for the D, D s , B, and 
B s systems. In Sec. II we establish the UMD principle 
first from experiment, then from a model calculation. We 
finally exhibit the UMD within a simple potential model. 
In Sec. Ill we use UMD to determine the parameters of 
the relativistic flux tube, a simple but fundamental model. 
These parameters are the Coulomb constant and the heavy 
quark masses m c and nib- Another application of UMD is 
given in Sec. IV where from the measured difference B s —B, 
a relationship between the constituent quark masses m s 
and m u ^d is established. This relationship is also shown 
to follow from relativistic kinematics alone. In Sec. V we 
use our results to predict a range of radially and orbitally 
excited HL mesons. 



II. UNIVERSAL LIGHT QUARK MASS 
DEPENDENCE 

The HL meson mass M , in the heavy quark limit, can be 
defined in terms of the excitation energy E and the heavy 
quark mass toq as 

M = rriQ + E . (3) 

As we will demonstrate, the meson mass M has universal 
mass dependence on both the heavy quark mass toq and 
the light quark mass to. Up to 1/toq corrections, we may 
expand E as 

E n>e = E n>e (0) + (3m 2 + . . . , (4) 
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where the coefficient /3 is independent of both the radial 
number n and the angular quantum number £. We ex- 
pect the expansion to have only even powers of m since 
m only appears quadratically in our model Hamiltonians. 
The energy differences between different HL excitations is 
then 



E 2 -E 1 = E 2 {0) - Ex(0) + /3(m 2 2 - to 2 ) + 



(5) 



The excitation energy differences of HL mesons with the 
same light flavor are independent of the quark mass. We 
offer three types of evidence for this (UMD) universality. 



A. Experimental Data 

We select any convenient P-wave and S-w&ve D type 
meson difference For example, 

Pi (2422 ± 2 MeV) - P(1864 ± 0.5 MeV) = 558 ± 2 MeV . 

We now compare this with the corresponding D s differ- 
ence, 

P s i(2535±0.5 MeV)-P s (1969±1.4 MeV) = 566±1.5 MeV. 

(7) 

If UMD is valid, the two differences should be identical. 
Experimentally they differ by 8 ± 3 MeV which is an ac- 
curacy of better than 2%. 

Other differences involving D* and D 2 give similar re- 
sults but with slightly larger errors. 



B. A dynamical model: the relativistic flux tube 

The Relativistic Flux Tube (RFT) (or QCD string) 
model with spinless quarks has been solved numerically for 
about a decade @. For a rigorous derivation and experi- 
mental motivation see We will not discuss the details 
of this model here except to emphasize that it is a very re- 
alistic model incorporating many of the features of QCD. 
In addition to the string confinement (with static tension 
a) we add a short range interaction U{r) = —k/r. In the 
heavy quark limit the heavy quark mass appears additively 
as in Eq. with no 1/toq corrections. The light quark 
constituent mass is not well known and we will only assume 
that it is in the range < to < 600 MeV. The parame- 
ters that appear in the RFT model are the string tension 
a, the Coulomb constant k, and the masses of the heavy 
quark, toq, and the light quark, to. Since the quarks are 
spinless in this model we will always compare our model 
predictions to the spin-averaged data. 

In this sub-section we are exploring the properties of the 
RFT model and not comparing to experimental data so the 
exact values of the parameters are not important. We use 
a = 0.18 GeV 2 and k = 0.5 which are in fact reasonable 
values, as we discuss in the next section. In Fig. 1 we 
plot the lowest S and P wave eigenvalues of the excitation 
energy E as a function of the light quark mass to. The 



important thing to notice is that the difference IP — IS 
is quite constant. This is exactly what is expected under 
UMD as in Eqs. 0) or @. 




FIG. 1: Light quark mass dependence of the lowest S and P- 
wave heavy-light excitation energies as predicted by the rel- 
ativistic flux tube model. We note that the IP — IS mass 
difference is nearly independent of light quark mass. 



C. A simple analytic result 

We show here that UMD is fundamentally a result of 
relativistic kinematics. Let us consider a simple time- 
component vector potential model with relativistic kine- 
matics, 



PTV = Eip, 



(8) 
(9) 



H = VP 2 + m 2 + U(r), 
where 

p 2 =p 2 r + L 2 /r 2 . (10) 

An expression of UMD is 

d 2 E/dL 2 dm 2 = 0. (11) 

We can demonstrate this to leading order with the 
Feynman-Hellmann theorem Q 



dE/dX = (dH/dX) , 



(12) 



where A is a system parameter. The desired quantity lllj l 
is then given by 

d 2 E/dL 2 dm 2 = (d 2 H/dL 2 dm 2 ) . (13) 

Expanding about L 2 = m 2 = we find from © that 

H=p r + (L 2 /r 2 +m 2 )/(2p r ) + ... (14) 
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which then yields to leading order 

d 2 E/dL 2 dm 2 = 0. (15) 

The above demonstration of light quark universality also 
holds for mesons with two light quark mesons. 

The next term in the expansion l|14|l is a cross term pro- 
portional to L 2 m 2 , which yields a non- vanishing second 
derivative l|15fl that violates UMD. A similar violation us- 
ing the RFT model can be seen in Fig. 1 and amounts to 
about 10 MeV for m increasing from 300 MeV to 500 MeV. 
This accounts for the small observed violation of 8±3 MeV 
noted in Eqs. © and 10. 



III. THE RFT PARAMETERS a, k, m c , m b 



by angular momentum multiplicity) states to be 

D sAS = ^D* s + ^D s = 2076 ± 1 MeV , (19) 

Ds.ip = ^D sfi+ + ±(Dl /2 + +D 3 J 2 + ) + ^D s , 2+ 

= 2515 ± 3 MeV . (20) 
We need the difference 

D s ,ip - D St is = 439 ± 4 MeV. (21) 

In Fig. 2 we show the RFT prediction for this difference as 
a function of the Coulomb constant k and we see that the 
correct value is, 

fc~ 0.52. (22) 



As we have noted the parameters entering the RFT 
model are the string tension a, the Coulomb constant k, 
and the two heavy quark masses m c and m^. The pre- 
dictions for excited states will be nearly independent of 
the light quark mass value but sensitive to the difference 
between m s and m u ^. 



A. String tension 

The universal Regge slope for both mesons and baryons 
is@, 

a' = 0.88 GeV 2 (16) 

For a relativistically rotating QCD string, the Nambu- 
Goto QCD string [|| and the RFT model predict the Regge 
slope to be 

a' = 1/2™, (17) 

which yields the string tension, 

a = 0.18 GeV 2 . (18) 

This value is quite consistent with that obtained from an 
analysis of heavy onia data alone and we will assume it 
in our subsequent work. 



B. Coulomb constant 

Without any assumptions about the quark masses, we 
can find the Coulomb constant using UMD as an organiz- 
ing principle. The idea is to compare the model predictions 
with an experimental IP — IS" HL mass difference. To do 
this we must know the masses of a pair of spin averaged 
states. Fortunately, we now have a complete set of states 
for the D s mesons due to recent B factory measurements 
and older data Q . We find the spin averaged (weighted 
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FIG. 2: The difference of the lowest S and P wave heavy-light 
meson masses in the RFT model as a function of the Coulomb 
constant k. The horizontal line is the experimental value (12 H 
determined from the D 3 states. 



C. Heavy quark masses 

The heavy quark masses do depend on the choice of light 
quark mass. In order to agree with the observed IS* state 
we must adjust rriQ as m is varied. The results for the 
charm and bottom quarks are shown in Figs. 3 and 4. 

IV. RELATION BETWEEN THE LIGHT QUARK 
MASSES 

For a given heavy quark, say the 6, the two ground state 
mesons B s and B are expected to differ in mass as given in 
Eq. © since /3 is not zero. Experimentally this difference 

is a 

B s - B = 91 ± 1 MeV . (23) 
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FIG. 3: The c quark mass required to yield the observed spin 
averaged Dis meson mass for a range of choices of light quark 
mass. 
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FIG. 4: The b quark mass required to yield the observed spin 
averaged Bis meson mass for a range of choices of light quark 
mass. 

This observation implies a functional relation between the 
strange and non-strange quark masses. Using the values 
for a and k in Eqs. I|18|) and l|22|) . we exhibit this relation 
for the RFT model in Fig. 5. We might note that the cor- 
responding charm difference is about 10 MeV larger and 
reflects a larger heavy quark kinetic energy (i.e., a 1/toq 
correction). Finally, we might comment that this relation 
between quark masses again arises primarily from relativis- 
tic kinematics. The relation 

yjp* + m\ - y/p 2 Q + ml d = 91 MeV (24) 

follows from the simple Hamiltonian ©. With the choice 
Pq = 0.4 GeV 2 as the average square momentum, the im- 
plicit relationship of Eq. Ij24(l parallels that of the more re- 
alistic RFT nicely and is depicted on Fig. 5 by the dashed 
curve. We also include the solutions for the light quark 



constituent masses obtained from analyses of hyperon mag- 
netic moments. We note that these "hyperon" values of the 
light quark masses fall close to our curve. 
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FIG. 5: The relation between the strange and the non-strange 
light quark masses required to maintain the relation 123H . The 
"data" points corresponds to the well known and successful 
hyperon magnetic moment model where the quarks have Dirac 
moments inversely proportional to the constituent quark masses 
(dot) and (error bar). The solid curve is the prediction 
of the relativistic flux-tube (RFT) model. The dashed curve 
is computed using the pure relativistic kinematics of Eq. 12 It 
with pi = 0.4 GeV. 



V. PREDICTIONS FOR THE HEAVY-LIGHT 
SPECTRA 

Now that the parameters of the RFT model have been 
fixed, we can make a range of predictions. We note that we 
have required only the S'-wave spin averaged ground states 
for the charm and bottom states, and one spin averaged 
P-wave multiplet (in our case the D s ). The predictions are 
then unique and independent of specific choices of the light 
quark masses. In Figs. 6 to 9 we present our predictions 
for the D,D S ,B, and B s flavor families. In each case we 
predict up to five radial and five angular states. As we 
expect from UMD, the predictions are nearly unique. If 
the light quark mass is varied over a 200 MeV range, the 
predictions for the excited states vary by less than 10 MeV, 
which would be difficult to see on the figures. In Table 1 
we provide the numerical predictions for the spin-averaged 
states assuming m„ i( j = 300 MeV and m s = 500 MeV. 

VI. CONCLUSIONS 

We have approached the subject of heavy-light meson 
spectroscopy by introducing a new principle which we call 
"universal light quark mass dependence" (UMD) . The idea 
is that the energies of all orbitally and radially excited 
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FIG. 6: The D(cu, cd) spectroscopy. The large solid dot is input 
data and the hollow dots are predicted states. The predicted 
states are displayed numerically in Table 1. The boxes repre- 
sent measured (spin) states not used in the calculation. The 
Di(2420) and D 2 (2460) are well known but the D*'(2637) [HI 
should be verified. 



FIG. 8: The B(bu, fed) spectroscopy. The solid dot is input data 
and the hollow dots are predicted states. The predicted states 
are displayed numerically in Table 1. 
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FIG. 7: The D s (cs) spectroscopy. The large solid dots are input 
data and the hollow dots are predicted states. The predicted 
states are displayed numerically in Table 1. To illustrate the 
spin dependence we show the IS and IP spin states as small 
dots. 



FIG. 9: The B s (bs) spectroscopy. The solid dot is input data 
and the hollow dots are predicted states. The predicted states 
are displayed numerically in Table 1. 



states vary in the same way as the light quark mass is var- 
ied. This proposal is supported in Sec.[n]by experimental 
evidence, numerical calculations using a realistic theoret- 
ical model, and finally by analytic demonstration using 
a simple but relativistic potential model. This universal- 
ity observation makes the analysis of heavy-light mesons 
transparent and considerably simpler. We further note 
that the measured B s — B difference implies a functional 
relation between the strange and non-strange light quark 
masses which is consistent with the well known quark 
model of the hyperon magnetic moments P, ||| . From only 
the 5-wave states and one spin-averaged P-wave state, the 



D s , we can reliably predict the D, D s , B, and B s excited 
spectrum. It should be noted that our predictions are for 
the spin-averaged states and that we assumed that the 
heavy-light assumption is valid. There are some small dis- 
crepancies in the data when thought of in the heavy-light 
limit. For example D s — D is about 10 percent larger (10 
MeV) than the corresponding difference B s — B. This is 
probably due to 1 /rriQ corrections to heavy quark symme- 
try. Another topic for future investigation is the break- 
down of the heavy-light approximation for highly excited 
states. 
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TABLE I: Predicted heavy-light meson states in GeV. The flux 
tube parameters are a — 0.18 GeV -2 , k = 0.524, m Ui d = 300 
MeV, m s = 495 MeV, m c = 1330 MeV, m b = 4670 MeV. 
These states are illustrated in Figs. 6-9. Varying the light quark 
masses over a range of 200 MeV changes our predictions by less 
than 10 MeV. 
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